EPJ manuscript No. 

(will be inserted by the editor) 



A Bayesian analysis of the nucleon QCD sum rules 

Keisuke Ohtani^'', Philipp Gubler^ and Makoto Oka^ 

Department of Physics, H-27, Tokyo Institute of Technology, Meguro, Tokyo 152-8551, Japan 
Received: date / Revised version: date 



o 

(N 
> 

o 
O 

m 

D 

& 

> 

vn 

^' 

o 



X 



Abstract. QCD sum rules of the nucleon channel are reanalyzed, using the maximum entropy method 
(MEM). This new approach, based on the Bayesian probability theory, does not restrict the spectral 
function to the usual "pole + continuum"-form, allowing a more flexible investigation of the nucleon 
spectral function. Making use of this flexibility, we are able to investigate the spectral functions of various 
interpolating fields, finding that the nucleon ground state mainly couples to an operator containing a scalar 
diquark. Moreover, we formulate the Gaussian sum rule for the nucleon channel and find that it is more 
suitable for the MEM analysis to extract the nucleon pole in the region of its experimental value, while 
the Borel sum rule does not contain enough information to clearly separate the nucleon pole from the 
continuum. 
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1 Introduction 

The use of QCD sum rules in studies investigating the 
properties of baryons has already a long history. Since the 
seminal papers of Shifman et al. [T], in which the QCD 
sum rule method was developed and the subsequent first 
application to baryonic channels by loffe [5] , the sum rules 
of the nucleon have been continuously improved by includ- 
ing higher orders in the perturbative Wilson coefficients 
[2113,^5,6,7,8, or non-perturbative power corrections [51131 
[TU] . This even lead to attempts to determine the mass dif- 
ference between the neutron and the proton which 
certainly is a very difficult task because of the smallncss 
of this difference compared to hadronic scales. Another 
important development was initiated in a paper by Lcin- 
weber |12| . in which, among other technical points, the 
choice of the interpolating field made by loffe was crit- 
icized and a new statistical method for the analysis of 
QCD sum rules was introduced. Furthermore, QCD sum 
rules also have been used to investigate the nucleon prop- 
erties in nuclear matter [T^[T^[T5] or at finite temperature 

However, in all these studies it was necessary to model 
the spectral function according to some specific functional 
form, the "pole -I- continuum" -ansatz being the most pop- 
ular one. Such a procedure inevitably incorporates strong 
assumptions on the spectral function into the analysis. 
This strategy works well when the actual spectral function 
has some resemblance with the chosen model, but will of 
course fail if this is not the case. For instance, as is known 



from studies using both QCD sum rules and lattice QCD, 
certain linear combinations of interpolating fields, which 
in principle carry the quantum numbers of the nucleon, 
couple to the nucleon ground state only very weakly. A 
QCD sum rule analysis of such interpolating fields, which 
uses the "pole + continuum" ansatz, can only lead to am- 
biguous results. The problems become even more severe 
in studies of the spectral function at finite temperature or 
density, as the validity of the phenomenological "pole -I- 
continuum" ansatz for the spectral function in such an en- 
vironment becomes less certain, rendering it more difRcult 
to make educated guesses about its actual form. 

To avoid the above mentioned problems, and to re- 
duce the assumptions that have to be made on the func- 
tional form of the spectral function, we have developed a 
novel method of analyzing QCD sum rules, which employs 
the maximum entropy method (MEM) |17| . The method, 
which is based on Bayesian probability theory, was shown 
to give reasonable results in the analysis of the sum rule 
of the p-meson channel. Moreover, recently this approach 
was successfully applied to the charmonium sum rules at 
both zero and finite temperatures |18j . 

In this paper, it is our main purpose to examine if and 
to what extent the MEM analysis can be applied to the 
sum rule of the nucleon. Throughout our investigations, 
we found that the MEM analysis of the Borel sum rule 



Gope(M) - 



dte 



'Pit), 
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in fact fails to satisfactorily extract the nucleon spectral 
function in the ground state region. As we will discuss 
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later in detail, this failure is mainly caused by the large 
contribution of the continuum to the OPE side of Eq.([T]), 
which strongly deteriorates the contribution of the nucleon 
pole. 



On the other hand, the Gaussian sum rule [TOIIM] 



Gope(s,t) 



dte' 



(2) 



which, for the nucleon, is formulated for the first time in 
this paper, turns out to give better results and allows us 
to resolve the nucleon pole from the continuum. There are 
essentially two reasons for the superiority of the sum rule 
of Eq.([2]). First of all, the kernel of Eg.®, a Gaussian cen- 
tered at s with a width of •\/2t, collects more information 
on the spectral function p(t) than the one in Eq.([T|) when 
the integration over t is carried out. This is especially true 
for small values of the width r, which we, however, cannot 
take arbitrarily small because below a certain threshold, 
the convergence of the operator product expansion (OPE) 
becomes poor. A similar situation occurs in the Borel sum 
rules, where the Borel mass is restricted from below due 
to the OPE convergence. The second reason also can be 
related to the kernel of Eq. ([5]) , containing two parameters 
s and T, which can be freely chosen as long as the OPE 
converges. This freedom allows us to vary two parameters 
at the same time in the MEM analysis of [T7], leading 
to reasonable results for the extracted spectral function. 
Similar experiences also have been made in nuclear struc- 
ture studies, where the Lorentz kernel has proven to be 
useful [null]. 

Furthermore, using the MEM analysis of the Gaussian 
sum rules, we are able to extract the spectral function 
not only in cases where the interpolating field strongly 
couples to the nucleon ground state and thus the "pole 
-|- continuum" ansatz should be valid, but also in cases 
where only higher energy states contribute to the sum 
rules and hence the conventional analysis most likely fails 
to give meaningful results. This advantage is especially 
useful for examining which kind of interpolating field cou- 
ples strongly to the nucleon ground state and is thus a 
suitable interpolating field for the analysis, a question with 
a long and controversial history in QCD sum rules studies 
[ gmm igS l lMllIl]. Our MEM analysis of the general op- 
erator given in Eq.® strongly suggests that the nucleon 
ground state only couples to f?i(x) (/3 = 0) and not to 
772 (x) (/3 = cxo)(see Eqs.Q and (gj) in Sec. [2]). In addition, 
we have obtained some hints of excited states coupling to 
772(2;). These issues will be discussed in detail in Section 

The paper is organized as follows. In Sec. [21 the details 
of the Borel and Gaussian sum rule for the nucleon are 
explained. Next, the maximum entropy method (MEM) 
is elucidated in Sec. [31 after which in Sec. [H the results of 
the analysis using the Borel sum rule in combination with 
MEM are presented. In Sec. [SI the results of the analysis 
for the Gaussian sum rule are outlined, and in Sec. 15.31 the 
differences of the obtained spectral functions depending on 
the choice of the interpolating fields are discussed. Finally, 
the summary and conclusions are given in Sec. [21 



2 QCD sum rules for the nucleon 

The method of QCD sum rules is used to carry out the 
analysis of the spectral function as follows. First, we choose 
an interpolating field which has the quantum numbers of 
the nucleon, then its correlation function is calculated in 
the deep Euclidean 4-momentum region. Alternatively, the 
same correlation function at the physical 4-momentum re- 
gion is expressed by the spectral function of the nucleon 
channel. The sum rules can then be constructed by equat- 
ing the two expressions using a dispersion relation. 

For the nucleon, there are two independent local inter- 
polating operators. 



mix) 
mix) 



'{u^''{x)C-i5d\x))u''{x) 



^abc ( ^,Ta 



{u"'{x)Cd\x))-i^u''{x). 



(3) 
(4) 

Here, abc are color indices, C is the charge conjugation 
matrix and T stands for the transposition operation. The 
spinor indices are omitted for simplicity. A general inter- 
polating operator can thus be expressed as 



ri{x) = 771(2;) + I3r]2ix), 



(5) 



where /3 is a real parameter. Here, the case of /3 = — 1 
is identified as the so-called "loffe current" [2], which is 
often used in QCD sum rule studies of the nucleon. 

Using this interpolating operator, we define the corre- 
lation function as 



n{q) = t / e'9^(0|T[77(2;)^(0)]|0)d^2; 



(6) 



^n,iq^) + n2iq^). 



The imaginary part of ni{q^) satisfies the positivity con- 
dition, lmlli{q^) > 0, while ImiT2(g^) is not necessar- 
ily positive due to contributions of negative parity states. 
The positivity condition is, however, essential for the ap- 
plication of the MEM method and we thus consider only 
i7i(g^) in the following. To make use of the information 
contained in Im772((7^), one would have to analyze the par- 
ity projected sum rules [32], which we plan to investigate 
in the future. In the deep Euclidean region (— — )■ ao), 
77i(g^) can be calculated by using the operator product 
expansion (OPE). Including operators up to dimension 8, 
we get 



niiq') 



5 + 2/3 + 5/ 3^ 
128(27r)4 



q'H~Q^) 



5 + 2/^ + 5/32 a 



256(27r)2 ' TT 
7-2/3-5/32 1 



^G') ln(-q2 



(7) 



24 q 

13-2^- 11;32 



96 



qq){qga ■ Gq) — 



noPEiq^ 



To obtain Eq. ([7]) , several approximations have been imple- 
mented. Firstly, only the lowest order in as is taken into 
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account. The validity of this approximation is not obvious, 
because it is known that the first order corrections are 
significant and lead to a considerable increase of the con- 
tinuum contribution [T2]. Nevertheless, our main goal of 
this paper is to examine whether the MEM analysis can 
be applied to the nucleon sum rule or not, and we thus 
ignore the as corrections here. For a more quantitative 
future analysis, the higher order corrections should cer- 
tainly be taken into account. The second approximation 
arises from the use of the vacuum saturation, by which 
(qqqq) and {qqqga ■ Gq) can be formally reduced to (qq)'^ 
and {qq)(qga ■ Gq), respectively. Although this approxi- 
mation can be justified in the large Nc limit [1], it is not 
clear to what extent it is trustable at iVc = 3. Nonethe- 
less, for the present qualitative analysis, we will assume 
this approximation to be valid. 

As already mentioned, ni{q^) can also be expressed in 
terms of the physical spectral function using the dispersion 
relation: 



Hliq') 



1 



Imi7i(<) 

t-g2 



dt 



Pit) 
t-q^ 



dt 



(8) 



^n,{q') 



where the definition ImTJi [t) = TTp{t) is used for the spec- 
tral function. Our goal is now to extract p(t) from the sum 
rule obtained by equating Eq.Q and Eq.(|8l). It should 
be noted here that subtractions are necessary in order 
to make the integral of Eq.® convergent. In the case of 
the nucleon, the subtraction terms are UiiO) + II[{Q)q^ + 

^^2*'°^ g^, which will disappear after transforming Eq.® 
into the Borel or Gaussian sum rules. How this is done 
will be explained in the following subsections. 



2.1 Borel sum rule 

In the case of the Borel sum rule, we transform IIi{q^ 
using the Borel transformation B^i, defined below: 



Bm = lim 



-,2,r.^^ _ 1)! 



(9) 



Applying Bm to Eq.(I7]), we get the following expression 
for Gope(M) = BmUoMi^)- 



64(27r)4 

5 + 2/3 + 513^ , a 



—G) 
256(27r)2 ' TT ' 

7-2^-5^ 2 1 
24 JP 

13-2/3- 11/32 



(10) 



96 



{qq){q9(^ ■Gq)j^. 



Meanwhile, applying the Borel transformation to Eq.®, 
we obtain for Gp{M) = B^npiq^): 



Gp{M) = 



A/2 



e ujp{Lj)diu. 



(11) 



Here , i = u;2 was used. Note that at this point, the sub- 
traction terms have been eliminated and the integral of 
the right hand side converges. This leads us to the final 
form of the Borel sum rule, 



Gope(A/) = Gp(A/) - — / e-^ujpiLo)dLu. (12) 



2.2 Gaussian sum rule 

The Gaussian sum rule, first introduced in |19j . exhibits 
another way of improving Eq.®. Based on the idea of 
local duality, it provides a formulation for the convolution 
of the spectral function with a Gaussian kernel. As this 
sort of sum rule is not often discussed in the literature and 
the specific case of the nucleon has to our knowledge not 
even been formulated, we will explain each step in some 
detail, following closely the formulation given in ^U\ . 

Going to the complex plane of q^ and taking the differ- 
ence between 7Ti(g2 = s + iA) and 7Ti(s — iA) of Eq.®, 
where s and A arc real, we obtain 



ni{s + iA) ~ ni{s ~ iA) 
2iA 



Pit) 



„ it-s)^ + A^^'- (1^) 

At this stage, the integral above is not convergent and the 
subtraction terms are not yet fully eliminated. Applying 
the following Borel transform B^r' 

(-Zi2V» 



and using 



we get 



B4 



BatI 



lim 



i—Y' 

{n-l)rdA^' 



1 (t-.)" 



(t-s)2 + zi2J 4r 



e 4^ 



(14) 



(15) 



iA) - 7Ti(s - iA) 



2iA 



1 



dte 



-pit) 



(16) 



V47rr Jo 

= Gp(s,T). 

Here, the subtraction terms have disappeared and the in- 
tegral in the above equation is convergent. 

As a next step, we will now show that Ea. (fT6)) can also 
be rewritten by using the inverse Laplace transform L~^. 
Using 

1 



(t - .S)2 + Z\2 



= L[e 



-x{t- 



"'^'dx, 



the kernel e can be re-expressed by i ^: 

1 



-x{t-sy 



1 

2Td 



(t-s)2+Z\2J 



b-\-ioc 



b — ioo 



1 



(17) 



{t~s)^ + A 



:dA' 
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Fig. 1. (a) Contour of integration A + /2 in Ea. (|18|) . The branch cut of IIi{—t) is represented by the wavy line on the negative 
real axis, (b) Closed contour C{R) of Ea. (|19p . The outer lines /a, /4 and Is are a part of the circle of radius R centered at —s, 
while r'e denotes a circular line with radius e, centered at the origin. The branch cut of 77i(— t) is represented by the wavy line 
on the negative real axis and the contour Fc is the straight line of the contour above and below the branch cut. The contours 
Fi and F2 are parts of Fi and F2, respectively. 



Setting X = the left-hand side of Eq. (IT51) ean thus be 
rewritten as 

1 _j^ ni{s + iA) - ni{s - iA) 



'47rT 



1 



1 



/Ittt 27ri Jf._ 



2iA ' 
ni{s + iA) - ni{s - iA) 
2iA 



]e^dA^ 



Then, replacing s + lA by —t in the first and s — iA by 
—t in the second term, we get 



Gp(s,T) 



1 



1 



/Attt 2-771 



77i(-t)exp(- ^ 7 )dt 
At 

(18) 

where the contour Fi + I2 is shown in Fig. [1^. 

Next, to obtain a sum rule that is practically usable, 
we consider the contour C{R) shown in Fig. [TJd. Taking 
the limit i? — > 00 and e — > 0, we are lead to the equation 
given below: 



Gp(s,T) 



= — lim 



1 



(t + s)"^ 



4t 



(19) 



Substituting the right-hand side of Eq.© into Ea. lfTO)) and 
examining the various terms, we see that the perturbative 
and dimension four terms only give contributions on the 
contour Fc- Meanwhile, the dimension six and eight terms 
do not have a branch discontinuity, but a pole at t = 
and therefore only contribute on F^. Using 



f exp(- 

J ~oo 



At 



)dt = V^[l+erf(^)] (20) 



where the error function erf (a;) is defined as 

2 



erf(a;) 



dx 



(21) 



and 



lim ■ 



we obtain 



Gope(s, t) 



±. ( lex {Jl±I^)dt 
27ri Jr^ F^^^^^ At ' 

i'-^'o {n - l)\ dt^-^^"^^^ At ' 
(n = 1,2,...), 

5 + 2/3 5/32 



(22) 



128(27r)4 

1,2 



xKl + ^)[l + erf(^)] + f^exp(-i-)] 

5 + 2/3 + 5^2 a g 
+ 512(2.)^ ^V^^t^+'-^^V^)^ 

= IS cxp(- — 

TTT 48 At 

1 13-2/3-11/32 



384 



{qq){q9o- ■ Gq)-exp{-—). 

T At 

(23) 

This then finally leads to the following form of the Gaus- 
sian sum rule, from which information of the spectral func- 
tion p{io) can be extracted: 

2 f°° {'^^-sf 
Gope(s, t) = Gp{s, t) = — / dw ■ uje ~- p(uj). 

vAttt Jo 

(24) 



'47rT Jo 

Here, we have again set t ^ lo'^ . 
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3 Maximum Entropy Method 

In this section, we will briefly explain the maximum en- 
tropy method (MEM). The general sum rule that we aim 
to analyze with this method is 



GoPE(a;) = Gp{x) 



dujK{x, aj)p(aj), 



(25) 



where K{x, to) stands for the kernel of either the Borel or 
the Gaussian sum rule, and x denotes the Borcl mass AI in 
the Borel sum rule or (s, r) in the Gaussian sum rule, re- 
spectively. Using Eg. ipsj) . we extract p{uj) from Gope(2^)- 
However, this is nontrivial, because GoPE(a;) contains a 
certain error a{x) arising from the uncertainty of the con- 
densate parameters and furthermore can only be used in 
the region where the OPE converges. For a specific value 
of X, Gp{x) is only sensitive to a limited region of /c(w), 
making it hard to obtain p{uj) in a wide range of lo. Be- 
cause of all these problems, rigorously solving Ea. ([?51) is 
in fact an ill-posed problem. Neverthless, the MEM tech- 
nique enables us at least to statistically determine the 
most probable form of p(a;). Using MEM has moreover 
the advantage that the solution of p{ijj) is unique if it ex- 
ists [57] ■ Therefore, different from the usual fitting, the 
problem of local minima does not occur. 

Let us now discuss the connection between MEM and 
Bayes' theorem. First, we define H to denote the prior 
knowledge on p{uj) such as positivity and asymptotic val- 
ues. Then. P[p\GH], the conditional probability of p given 
GoPE and H is rewritten using Bayes' theorem as 



P[p\GH] = 



P[G\pH]P[p\H] 
P[G\H] 



(26) 



Here, P[p|i7] is the so-called prior probability, and P[G\pH] 
stands for the likelihood function. The most probable form 
of p{(jj) is obtained, by maximizing P[p\GH] of Ea. ((26|) . 
For carrying out this task, one needs the concrete func- 
tional forms of P[/o|i/] and P[G\pH]. which can be ex- 
pressed as pniST] 



P[G\pH] = e 
P[p\H] 



(27) 



where a is a real positive number and the functionals L[p\ 
and S[p] are defined as 



L[p]^ 



1 



dx' 



and 



a^x) 



(28) 



S[p] = / duj[p{uj) - m{cu) - pioj) log{^)], (29) 
Jo 'm-l^) 

where, the error (t{x) in Eg. ([55)1 is evaluated by the sta- 
tistical method explained in [12]. S[p] is also known as the 
Shannon-Jaynes entropy and the positive function m(a;) 
appearing in this expression is called the default model. 



Since P[G|i/] does not depend on p{oj), we drop it as 
it corresponds only to a normalization constant. Hence, 
P[p\GH] is given as 

P[p\GH] (X P[G\pH]P[p\H] 

- pQ[p] 



where, 



Q[p] = aS[p] - L[p], 



(30) 



which means that to get the most probable p{uj), one has 
to solve the numerical problem of obtaining the form of 
p{uj) that maximizes Q[p], for which we use the Bryan 
algorithm |28| . Once a Pa{^) that maximizes Q[p] for a 
fixed value of a is found, this parameter is integrated out 
by averaging p over a range of values of a and assum- 
ing that P[p\GH] is sharply peaked around its maximum 
PW\GH]: 



PoM = J [dp] J dap{uj)P[p\GH]P[a\GH] 



(31) 



dapa{oj)P[a\GH] 



The integrated range of a is P[pc I G77] > 0.1xP[p„„,JG7J] 
where ofmax is the value of a which maximizes P[pa\GH]. 
By using Bayes' theorem, P[a\GH] is expressed as (for 
details see [?7]): 

P[a\GH] cx P[a\H] cxp(i V log + Q[p„]), 

2 ^ a + Xk 



where are the eigenvalues of the matrix. 



dpidpj 



and Pi denotes p(wi) dcu. 

In the numerical analysis, Eas. (|25l28l29p are discretized 
as follows: 



G{xi) ^ ^ duK{xi,ujj)p{ujj), 



(32) 



L{P) 



2{x 



max •^min J 



[GopE(a;i) - Gp{xi 



a'^{xi) 



(33) 



Pj^j) - 
m{ujj)' 



(34) 

In Eas. ([5^ - we take 100 data points in the analyzed 
region of the variable x (A'a;=100) and 600 data points for 
the variable lo (A'„=600). Furthermore, we adjust cjuiim 
duj and dx to 0[GeV], 6[GeV], (wmax - ^mm)/N^ 
and (xmax - Xniin)/Na:, respectively. 

Finally, the method of evaluating the error of the ob- 
tained spectral function is introduced below. We estimate 
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the averaged error of Pout{^) over a certain region [wi, a;2]. 
{{Sp a)'^)uji,uj2i the dispersion of Pa{(^) in [wi,a;2] is defined 
as 

= 7 f[dp] f ' dujdLj'Sp^{Lu)dpa.{uj')P[p\GB] 



(W2 - t^l)^ ^1 



dujduj'{ 



6p{u!)Sp{uj') ' 



lP=Pa ) 



where Spa{uj) = p(ui) — pa{ui) and the gaussian approxima- 
tion of P[p|GH] around p^ has been used. After averaging 
over a in the same way as in Ea. pip . we arrive at the final 
value of 5 Pont ■ 



(<5pout)wi, 



day^((<5p„)2)^^,^,P[a|GH]. (35) 



The errors of the spectral functions shown in the figures of 
this paper were all estimated by this method. These errors 
will be indicated by three horizontal lines, whose lengths 
stand for the region [wi , 102] , over which the average is 
taken, while their heights correspond to (pout) + i^Pout), 
(Pout), (pout) — {S Pout) , respectively. For more details about 
MEM, we refer the reader to [7fll29] . For an application 
of this method to the nucleon channel in the framework of 
lattice QCD, see [3D] and for a discussion of specific issues 
related to QCD sum rules, consult [T7]. 



> 

o 

O 
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1 1.1 1.2 1.3 

Borel mass M[GeV] 



1.4 



Fig. 2. Comparison of the mock data and the OPE data at 
13 = -1. The solid line shows Gope(M) of Eq.dTOj, with the 
respective error indicated by the gray region. The dashed line 
denotes Gmock(Af) for which the mock spectral function of 
Ea. p6p was used. 







{qga ■ Gq) / (qq) 


-(0.24 ± 0.01)^ GeV^ 


0.012 ± 0.0036 GeV* 


0.8±0.2GeV'' 



Table 1. Values of the parameters appearing in the OPE, 
taken from 1311. 



where we use the following values for the various param- 
eters: 



4 Analysis using the Borel sum rule 

In this section, we will analyze the nucleon spectral func- 
tion for the Borel sum rule. It is easily understood from 
dimensional considerations that unlike in the meson case, 
the contribution of the continuum states to the baryon 
spectral function is proportional to ui* and thus strongly 
enhanced. As was done in similar studies using MEM and 
lattice QCD, we will therefore analyze p(w)/a;^ instead 
of p{uj) and hence from now on denote p{u})/uj^ as p(w), 
leading to the equations below: 

/•OO 

Gope(A/) = Gp(A/) = / du;K{M,uj)p{uj), 
Jo 



2^5 



4.1 Analysis using mock data 

In order to check the effectiveness of MEM to extract the 
spectral function of the nucleon, we first carry out an anal- 
ysis using mock data. The employed mock spectral func- 
tion is given below: 



ock 



2M^ 



S{u}-Mn) 



5 + 2P + 



1 



128(27r) 



Mn = 940 MeV, uq - 
6 = 0.05 GeV, 



1.3 GeV 
0.19 



(2^) 



GeV' 



(37) 



Here, Pmock('^) is constructed to have a narrow ground 
state pole and a continuum, which approaches the pertur- 
bative value at high energy. Defining now 



Gmock(Ajf) 



dujK{M,uj)p,^ock{(^), 



(38) 



we apply the MEM procedure to 



G 



mock 



(M) ^ Gp{M) 



dujK{M,uj)p{u). (39) 



The residue of the nucleon pole A was fitted so that Ginock(Af) 
matches Gope,{M) at /3 = — 1 in the analyzed Borel mass 
region. For comparison, GoPE(Af) and Gp{M) are shown 
in Fig.[2j The values of the parameters appearing in Eq.© 
are shown in Table[TJ As can be observed in Fig. [31 Gopb{M) 
and Gp{M) are consistent within the range of the er- 
ror, which shows that the "pole -I- continuum" ansatz de- 
scribes the OPE data well. However, this does not neces- 
sarily imply that the integral of Eg. ([5^ can be reliably 
inverted and that valid information on the nucleon pole 
can be extracted. To investigate to what extent this is 
possible, we now analyze Gniock(A^)- For a realistic anal- 
ysis, we here employ the error obtained from Gope(.^^), 



(36) crmock(M) = croPE(Af). 
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Fig. 3. Spectral functions extracted from Gmock(Af) using various default models at /3 = —1. The analyzed Borel mass region 
is 0.91 GeV < M < 0.97 GeV. The solid line shows p(cj), the dashed line depicts the input mock spectral function and the dash- 
dotted line stands for the default model. As we explain in Section 3, the three horizontal lines correspond to (pout) + (Spout), 
(pout), (Pout) — (5pout), respectively. Similar horizontal lines in following figures are understood to have the same meaning. 



To use MEM, we must at first fix the default model 
m{u!). A reasonable choice for 771(0;) should rcficct our 
prior knowledge on the spectral function such as the asymp- 
totic behavior at low or high energy. To test several possi- 
ble choices, we here introduce three types of default mod- 
els. The first one is a constant consistent only with the 
asymptotic behavior of the spectral function at low en- 
ergy, therefore lying close to 0: 



mflat(w) = 1.0 X 10" 



(40) 



The detailed value of msat is not so important, as long 
as it can be considered to be small enough compared to 
the asymptotic value at high energy. As we will discuss 
in more detail in the section dealing with the gaussian 
sum rule, we indeed have found that the position of the 
lowest lying pole and its residue of the obtained spectral 
functions depend on the value of Wfiat only weakly. The 
second default model is also a constant which now reflects 
the asymptotic behavior at high energy: 



5 + 2/3 + 5(3^ 
128(27r)4 



(41) 



The third one is a combination of the first two with the 
correct behavior at both high and low energy: 



TOhybr(w) = 



5 + 2/3 + 5/32 



1 



128(27r)4 
Wo = 1.8 GeV, (5 = 



1 _(_ g(wo-w)/<5 ' 

0.1 GeV. 



(42) 



In the following, the analysis is carried out at /3 = —1. 
In the case of the conventional method, which assumes the 
spectral function to have a specific functional form such as 
the "pole -|- continuum" ansatz, the analyzed Borel mass 
region is restricted so that the ratio of the highest dimen- 
sional term is less than 0.1 of the whole Gope(-^^) to have 



some confidence on the OPE convergence. The Borel mass 
region is further limited by the condition that the contri- 
bution from the continuum states in Gp{M) should be less 
than 0.5 to make sure that the lowest pole dominates the 
sum rule. The Borel mass region determined according to 
these restrictions becomes 0.91 GeV < M < 0.97 GeV. As 
pointed out in |12j , this region is very narrow and thus we 
expect that it will be very difficult to extract p{u;) in a 
wide range of uj with this small amount of available infor- 
mation. Although, when using MEM, we can in principle 
employ values of M above 0.97 GeV, we here first analyze 
the spectral function using the Borel mass region 0.91 GeV 
< M < 0.97 GeV. The resuhs are shown in Fig.E] 



It is clear from Fig. [3j that the obtained lowest peaks 
lie much above the input value of 940 MeV. Hence, as ex- 
pected, we cannot extract much information on the nu- 
cleon pole from this sum rule. Especially, in the case of the 
default model masym(w), the spectral function in the low 
energy region approaches TOasym(O) and we can only ob- 
serve a small lowest peak. A similar tendency was observed 
in the p-meson channel [T7]- We will therefore abandon 
this default model in the following. In the case of TOflat(w), 
although the high energy behavior wrongly approaches 
at high energy, the low energy behavior, which is the main 
focus of our interest, seems to be reasonable. On the other 
hand, using mhybr(w) leads to the correct behavior at both 
high and low energy. From these results, we can infer that 
the default model completely determines the asymptotic 
values of the spectral function. This behavior at high and 
low energy should therefore be considered to be an input 
in the current analysis. This can be understood from the 
properties of the kernel K{M,uj), leading to a function 
Gp{M), which is insensitive to the values of p(w) at large 
and small values of w. From the behavior of K{M,u!), we 
can also expect that increasing the upper boundary of M 
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Fig. 4. Spectral functions extracted from Gniock(Af) using var- 
ious default models at /3 = — 1. The analyzed Borel mass region 
is 0.91 GeV < M < 1.4 GeV. The solid hne shows p{lj), the 
dashed line depicts the input mock spectral function and the 
dash-dotted line stands for the default model. 
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Fig. 5. Spectral functions extracted from Gope(M) using var- 
ious default models at /3 = —1. The analyzed Borel mass re- 
gion is 0.91 GeV <M< 0.97 GeV. The solid line shows p(uj), 
the vertical dashed line indicates the position of the nucleon 
ground state and the dash-dotted line stands for the default 
model. 



allows the analysis to become sensitive to p{u>) at higher 
energy regions. 

For investigating this case, we analyze the spectral 
function under the condition of Mn^n = 0.91 GeV and 
-^max = 1.4 GeV. The results are shown in Fig. 01 When 
using mflat(w), the resulting spectral function at high en- 
ergy oscillates around the continuum value before approach- 
ing the default model. This is plausible, because the input 
spectral function has the "pole + continuum" structure, 
which the IVIEM analysis is trying to simulate with the 
given limited information of Gmock(-M). Nevertheless, the 
MEM procedure cannot reproduce the mass of the nu- 
cleon, whatever default model or Borel mass range is used. 



4.2 Analysis using OPE data 



Similar to the previous section, we now carry out the anal- 
ysis using the real OPE data, Gope{M), even though from 
our experience of the mock data analysis, we cannot ex- 
pect to obtain meaningful results. We analyze the spec- 
tral function by setting the analyzed Borel mass region to 
0.91 GeV < M < 0.97 GeV, as in the last section. The re- 
sults are shown in Fig. [SJ Using a wider Borel mass region 
leads to spectral functions very similar to the ones shown 
in Fig. m As in the mock data analysis, the MEM proce- 
dure does not succeed to reproduce the nucleon peak. The 
main reason for this failure can be traced back to the slow 
convergence of the OPE and to the large contribution of 
continuum to the sum rule. These factors severely reduce 
the information of the lowest nucleon pole that can be 
extracted from C?oPE(Af). 



5 Analysis using the Gaussian sum rule 

In case of the Gaussian sum rule, the analysis is carried 
out using the following equations: 

poo 

Gope(s, r) = Gp(s, t) = / dujK{s,T,uj)p{uj), 

Jo 

V47rr 

As before, we set /3 = — 1 and the results will be shown in 
terms of the dimcnsionless spectral function. 

5.1 Analysis using mock data 

As for the Borel sum rule, we use Eg. ((36)) as mock data and 
Eqs. (PDl - H^ for the default model. As a first step, we must 
determine the ranges of s, r used in the analysis. From 
the property of the kernel, we expect that G p at small 
values of t will be more sensitive to narrow structures 
such as the lowest peak, while Gp at larger values will to 
a large extent be fixed by the continuum. Hence, to extract 
as much information as possible from Ginock(s, t), wc use 
several value of t at the same time, which are t = 0.5, 1, 
1.5, 2 GeV''. From 

lim ^^e"*"'*-°'' ^ b(u? ~ s), (43) 
T^o ^47rr 

one understands that in the limit t — > 0, Gp{s, r) in prin- 
ciple should approach the spectral function as 

limGp(,s,r) = sV(v^)- (44) 

However, it is not possible to take this limit, because the 
OPE is an expansion in powers of 1/ \/t as can be seen in 
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Fig. 6. Three examples of the kernel K{s,t,ijj). The solid line 
shows K(s,T,uj) for s — — 0.79GeV^, the dashed line for s — 
GeV^ and the dash-dotted line for s = 0.94 GeV^ In all three 
cases, r is fixed to r = 0.5 GeV*. 



T [GcV**] 


0.5 


1 


1.5 


2 


Smin [GCV"] 


-0.79 


-1.84 


-2.96 


-4.13 


Srnax [GeV"] 


0.94 


0.02 


-0.97 


-2.0 



Table 2. The analyzed ranges of the variable s for the em- 
ployed values of t. 



Ea. (j23p . meaning that the convergence of the expansion 
worsens with lower values of t. 

Turning now to the parameter s, similar to the Borcl 
sum rule case, we determine its range from the conver- 
gence in Gope(s,t) and the contribution from the con- 
tinuum in Gp{s,T). The analyzed regions of s and r are 
shown in Table [5J Here, we denote the upper and lower 
boundaries of s at each r as Smax and Smin, respectively. 

Some further comments are in order here. Firstly, let 
us briefly explain the method of using the two variables of 
s and T at the same time, compared to the discussion of 
Section 131 where only one variable x was considered. Gen- 
eralizing the method from one variable to two is straight- 
forward, as one only needs to redefine the kernel and the 
likelihood function. Specifically, we use 



G(a;j) = ^ K{xi,LLij)p{ujj)duj 

CJj — CJmin 



2(^max •^min) _ 



J2 dx 



aHx,) 



r*^— Tmin Si—S 



^(-^max "^min) 

In the present analysis, we take 4 data points for t [N^ = 
4), and 25 data points for s [Ng = 25) at each and 
adjust Wmin and Wmax to and 6 GeV, respectively, s^^x 
and s^ijj at each are given in Tabled and ds^ is (s^^x^ 

Secondly, we comment on the ranges of the variable s, 
shown in Tabled] It is seen that to the most part, we use 
negative values for s, which may seem somewhat counter- 
intuitive as, naively, the kernel of Ea. (P5)) seems to have 
a peak at y/s, when considered as a function of uj. There- 
fore one would expect values of s around 1 GeV to give 
a Gp{s,T) that is most sensitive to the spectral function 
in the region of the nucleon pole. However, this is not the 
case because the kernel is distorted due to the w^-factor 
in front of the Gaussian of Eg. pS)) . For example, for con- 
structing a kernel with a maximum value at 1 GeV which 
is the energy region which we are mostly interested in, s 
becomes 1 — 5t and, in the case of t = 0.5, 1, 1.5, 2 GeV^, 
has a negative value. This is illustrated in Fig. [51 where 
one can see that only kernels with negative s have peaks 
around 1 GeV. From these arguments, we can understand 
that it is necessary and important to use negative values 
of s in the analysis of the Gaussian sum rules. 

Let us now discuss the obtained spectral functions, 
which are shown in Fig. [71 Their detailed numerical re- 
sults are given in Table [31 It is observed that, compared 
with the Borel sum rule, the reconstruction of the lowest 
lying peak has considerably improved. Nevertheless, it is 
seen in Table [31 that its position is shifted upwards about 
120 ~ 180 MeV, which should give an idea of the preci- 
sion attainable with this method. To evaluate the residue 
of the first peak, we first have to define the region ui over 
which the peak can be considered to be the dominant con- 
tribution. The left edge of the first peak is determined to 
be the point whose height is 1 /30 of the peak vertex and 
the right edge so that the center between the left and the 
right edge lies on the peak vertex. The residue is then ob- 
tained by simply integrating p{uj) over the peak region, 
indicated by the horizontal length of the error bars. Table 
[31 shows that the residue of the first peak gives about 80 
% of the input value of the mock data. Since the default 
model TOflat('^) leads to a peak position and residue closest 
to the input values we will use only TOfiat(w) as the default 
model in the following. 

Here, we comment on the width of the peaks appear- 
ing in the obtained spectral functions. One might wonder 
why these peaks have a finite width even though the input 
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Fig. 7. Spectral functions extracted from Gmock(s,r) using various default models at /3 = —1. The analyzed regions of s and r 
are shown in Table [J] The solid line shows p(cj), the dashed line depicts the input mock spectral function and the dash-dotted 
line stands for the default model. 





mock data 




mhybr 




position of first peak [MeV] 


940 


1060 ± 140 


1120 ± 100 


1080 ± 410 


position of second peak [MeV] 




1850 


1970 


1960 


residue of first peak [lO^'^'GcV] 


8.3 


7.0 


6.0 





Table 3. Position of the peaks and values of the corresponding residues, obtained from a mock data analysis with three default 
models at /? = —1. Here, the given values of the mock data residue is defined as \^/{2M%). The errors are estimated by 
evaluating the half width of the obtained first peak. 



mock data only contains a zero width pole. There are two 
important aspects related to this point. Firstly, one should 
note that the OPE side of the sum rule is rather insensi- 
tive to the value of the width, because this information is 
largely lost after the integration over the spectral function 
in Eqs. (IT|) or This is so as long as the parameters gov- 
erning the scales of these integrals (M or r^^^) are much 
larger than the width of the peak of interest. This point 
was discussed explicitly for the p- meson channel in |17j . 
Secondly, one should remember that MEM is a statistical 
method, that can only provide the the most probable form 
of the spectral function, given the information available. 
This most likely spectral function depends not only on the 
input data but also on their error. Generally, a larger error 
makes a peak broader and smaller in magnitude, therefore 
introducing an artificial width connected to the error of 
the OPE data. 

Next, to study the dependence of the results on the 
range of s, we fix s^ax to 1.0 GeV^ 1.5 GeV^ and 2.0 GeV^ 
for all values of r, and redo the analysis. The results are 
shown in Fig.[5]and Table|31 From Tabled one can observe 
that the position and residue of the first peak depends on 
Smax only weakly, while the position of the second peak is 
quite sensitive to the value of Smax- Furthermore, from the 
Fig. [51 it can be understood that the second peak rather 



represents the continuum, around which the MEM output 
oscillates. 



Finally, we discuss the dependence of the spectral func- 
tion on the parameters appearing in the default models. 
In the case of the flat default models as given in Egs. pn]) 
and (|¥T|) , we have carried out analyses using several height 
values (spanning over several orders of magnitude). The 
results are shown in Fig.[S] It is clear that the detailed val- 
ues have little influence on the position of the nucleon pole. 
Extracting the corrcspoding residues, it is understood that 
their values also depend on the height value only weakly. 
Specifically, the pole positions for the curves shown in Fig. 
[9] coincide within 1% and the residues within 6%. In the 
case of a hybrid default model as given in Ea. (|^ . we 
also carried out analyses using several parameter combi- 
nations, the result being that the positions of the obtained 
first peaks did depend on the parameter combinations only 
weakly. 

In all, we conclude that using the Gaussian sum rule, 
reconstruction of the lowest lying peak and the residue of 
mock data is largely improved compared to the Borel sum 
rule, and that the attainable precision with the present 
uncertainties of the condensates is about 20 %. 
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Fig. 8. Spectral functions extracted from Gmock(s,''") and mnati^) at /3 = —1. Various values of Smax are used as indicated in 
the figures. The values of Smin at different r are shown in Table [J] The solid line shows p(tj) and the dashed line stands for the 
input mock spectral function. 



Smax [GeV"] 


Smax — 1.0 




Smax — 2.0 


position of first peak [MeV] 


1090 ± 140 


1100 ± 140 


1100 ± 150 


position of second peak [MeV] 


1980 


2020 




residue of first peak [10~^GeV] 


6.7 


6.5 


6.3 



Table 4. Position of the peaks and values of the corresponding residues. Here, the default model mflat(i^) and various s- 
values were used. The errors are estimated by evaluating the half width of the obtained first peak. 



3 



2 









A 
i"l 
i,i 

n 
H 


default model = 1 .Ox 1 0"^ 

default model =1.0x10"^ 

default model =1.0x10 




/'i \ 

k i 
a !\ 

/// '>V_^ 


A 

— 









12 3 

0)[GeV] 



Fig. 9. The spectral functions ectracted from Gmock(s,r), 
using various height values for the fiat default model mflat((^). 
The concrete numerical values are given in the figure. 



5.2 Analysis using OPE data 

Next, we carry out the analysis using GopeCsjt). We ap- 
ply the MEM in the same regions of s and r used in the 
mock data analysis. The results are shown in Fig.[TUJ while 
the numerical details are given in Table [51 The behavior 



of the results is similar to those of the mock data analysis. 
We observe that the positions of the lowest peak lie quite 
close to the experimental value. Besides a clearly resolved 
first peak, the spectral function exhibits one or two higher 
peaks which oscillate around the asymptotic high energy 
limit. In principle, this part could also contain nucleon res- 
onances with both positive and negative parity. With the 
present resolution achievable with the MEM technique, it 
however seems to be rather difficult to resolve these res- 
onances from the continuum. Using the parity projected 
sum rules [22], which can also be analyzed by the MEM 
approach, might improve the situation. 

5.3 Investigation of the /3 dependence 

The coupling strengths of the nucleon ground state and its 
excited states depend on the choice of the interpolating op- 
erator, i.e., on the value of (3. To investigate the nature of 
this dependence, we let /3 vary as —1.5, —0.5, 0.0, 0.5, 1.0, 
1.5, oo and extract the corresponding spectral functions. 
For r, we have chosen 0.5, 1.0, 1.5, 2.0 GeV* as before and 
Sniin is determined from the OPE convergence condition at 
each /3. In the case of /3 = 1, Smin is determined from the 
dimension 4 term because the higher dimensional terms 
vanish for this choice. To obtain information not only on 
the spectral function around 1 GeV, but also in the region 
of possible excited states, s,nax is not determined from the 
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Fig. 10. Spectral functions extracted from the MEM analysis using Gope(s,t) and mflat('^) at /3 = —1. The solid lines show 
p{oj) and the dash-dotted lines stands for the asymptotic value of the spectral function at high energy. (Note that p{uj) is 
normalized by ^.) In the top figure (cont), the values of s and r given in Tableware used, while for the other figures the values 
of Smax are fixed as indicated. 



Smax [GeV"] 


Smax = (cont) 




Smax — 1 - 5 


Smax — 2.0 


position of first peak [MeV] 


990 ± 130 


1040 ± 140 


1060 ± 150 


1050 ± 150 


position of second peak [MeV] 


1840 


1920 




1800 


residue of first peak [10~^GeV] 


7.9 


7.0 


6.6 


6.3 



Table 5. Positions of the first and second peaks and the corresponding residues of the first peak obtained from the analyses 
using OPE data at /3 = — 1 and the default model mflat(t»j). The errors are estimated by evaluating the half width of the obtained 
first peak. 



/3 


-1.5 


-1 


-0.5 





0.5 


1 


1.5 


oo 


r =0.5 


-0.40 


-0.79 


-0.39 


-0.35 


-0.43 


0.92 


-0.23 


1.04 


r =1.0 


-1.21 


-1.84 


-0.96 


-0.87 


-1.10 


0.12 


-0.93 


-0.36 


T =1.5 


-2.12 


-2.96 


-1.63 


-1.47 


-1.87 


-0.76 


-1.74 


-0.88 


T =2.0 


-3.07 


-4.13 


-2.34 


-2.13 


-2.68 


-1.66 


-2.60 


-1.45 



Table 6. Values of Smin [GeV^] at each /3 and t [GeV*]. s^ax is fixed to 2.0 GeV^. 



dominance of the lowest lying state in Gp(s, r) but is fixed 
to 2.0 GeV^. The explicit values of Snun are given in Table 
[HI Note that /3 = oo should not be taken literally, but just 
means that we use the correlator of only 772(2;) (which is 
the coefficient of in Eq.(I7])) for the analysis. The result- 
ing spectral functions and their numerical properties are 
shown in Fig. [TT] and TablejT] For a better comparison, we 
have normalized the spectral functions in Fig.[TT]by divid- 
ing them by the factor TV (^) = {5 + 2 13 + 5 fi^)/ {128 {2n)'^), 
so that the continuum approaches unity in the high energy 
limit. 

Considering first the lowest peaks at w ^ 1 GeV, we sec 
that a peak is clearly resolved for /3 = —1.0, —0.5, 0.0 and 
0.5, located at about the same position and with similar 
residue values. For the other values oi 13, no prominent 



peak is observed. These results are, however, obtained by 
using the Smin values determined separately for each /3 
and one could suspect that the choice of Smin affects the 
properties of the first peak. To study this dependence and 
to get an idea of the stability of our results, we redo the 
analysis, now using values of Smin that are determined via 
a /3 independent criterion as advocated in [T^[55) . Taking 
a closer look at Table HI it is observed that in the region 
between f3 = —1.0 and 0.5, where the ground state peak 
can be extracted, the OPE convergence is worst for f3 = 
0.0, giving the largest value of Smin- This implies that it is 
most reasonable to fix Smin to the values of /3 = 0.0, so that 
no values of s are used, where the convergence of the OPE 
might be questionable. Using this criterion, we repeat the 
analysis at (3 = —1, —0.5, 0.5. As before, we set Smax 
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Fig. 11. Spectral functions obtained from the MEM analyses of OPE data for various values of /3, normalized by the factor 
N{l3). The values Smin are given in Table [6] and Smax is fixed to 2.0 GeV^. 





-1.5 


-1 


-0.5 


0.0 


0.5 


1 


1.5 


oo 


position of first peak [MeV] 


1210 


1050 


1120 


1130 


1120 


1300 


1470 


1630 


residue of first peak [10~^GeV] 




6.3 


5.0 


4.8 


5.3 









Table 7. Positions of the first peaks of the spectral functions shown in Fig. 1111 For the cases in which the first peak can be 
interpreted as the nucleon ground state, the corresponding residue is given as well. 



to 2.0 GeV^. The results are shown and compared to the 
previous ones in Fig. 1121 The numerical details are given in 
Table El We see from these results that some details of the 
output spectral functions (especially at /3 — —1) change 
when the new criterion is used. The qualitative structure 
of a clear lowest peak with a continuum structure at a 
somewhat higher energy is however not altered. 

We also observe from Table [8] that the position and 
residue are almost independent of /?. This fact becomes 
even more explicit when we plot the non-normalized spec- 
tral functions for several values of f3 around /3 = as 
shown in Fig. [T31 This figure clearly illustrates that, in 
contrast to the continuum spectra, the property of the 
lowest peak is essentially /?- independent. 

All these results can be naturally interpreted by pre- 
suming that the lowest lying pole couples strongly to r]i, 
and only very weakly to 772- Therefore, the pole appears 
to be almost independent of /3 as long as |/?| is small, 
but disappears when becomes large and the contribu- 
tion of 772 dominates the spectral function. The nucleon 
ground state pole can be resolved as long as its strength 
is large enough compared to the continuum contribution. 



At about > 1.5, however, this continuum contribution 
gets too large to extract information on the lowest peak. 
These conclusions agree with the findings of lattice QCD 
^3,34,35] and some earlier suggestions of QCD sum rule 
studies ^4], in which however the position of the nucleon 
ground state peak was used as an input. 

Furthermore, by looking at Table E] and Fig. [T21 we 
can also understand why it is the lofFe current {f3 = —1.0), 
rather than rji, that seems to be most suitable for QCD 
sum rule studies. This comes from the fact that the OPE 
convergence is considerably faster for the loffe compared 
to rji , which allows the analysis to become more sensitive 
to the lower energy region of the spectral function, there- 
fore providing more information on the nucleon peak. 

As a last point, let us also consider the continuous 
structure above the nucleon ground state pole and possi- 
ble excited states appearing there. It can be observed in 
Fig. [TT] that for most values of /?, the spectral function 
just oscillates around the high energy limit, similar to the 
results obtained from the mock data analysis shown in 
Fig. IHl where we have only included the continuum and 
no resonances into the input spectral function. On the 
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Fig. 12. Spectral functions obtained from a MEM analysis of OPE data with mfiat(i^) at /3 = —1, —0.5 and 0.5, normalized 
by N(j3). The solid lines show the normalized output spectral function for which a /3 independent "refined" criterion was 

used to determine Smin. The details are explained in the text. For comparison, the spectral functions obtained with the previous 
/3-dependent criterion are indicated by the dashed line. 



/3 


-1 


-0.5 





0.5 


position of first peak [MeV] 


1120 


1120 


1130 


1140 


residue of first peak [10^'' Go V] 


5.3 


4.9 


4.8 


5.1 



Table 8. Positions of the first peaks and the corresponding 
residues of the spectral functions shown in Fig. 1121 where a 
/3-independent criterion was used to fix Smin. For comparison, 
the result for fi—Q is shown as well. 
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Fig. 13. Non-normalized spectral functions obtained from a 
MEM analysis using OPE data and mflat(t<j) at 13 — —1.0, —0.5, 
0, 0.5. The employed Smin is fixed to the one at /3 = and Smax 
to 2.0 GeV^ 



other hand, for /3 = 1.5 and especially /3 = oo, the be- 
havior is somewhat different, showing a quite clear first 
peak at about 1.6 GcV. This peak could correspond to 
one of the nucleon resonances, N(1440), N(1535), N(1650) 



or a combination of them. Results of lattice QCD stud- 
ies [M1I55] indicate that 772 couples to the negative parity 
states N(1535) and N(1650), while the first excited posi- 
tive parity state lies considerably above the Roper reso- 
nance N(1440). Even though the systematic uncertainties 
of these calculations are probably not yet completely un- 
der control, these findings suggest that the lowest peak 
seen at /3 = 00 corresponds rather to N(1535) or N(1650) 
and not to N(1440). It, however, seems to be difficult at 
the present stage to make any conclusive statements on 
the nature of this peak. An analysis of the parity pro- 
jected sum rules will hopefully help to clarify this issue. 



6 Summary and Conclusion 

In this study, we have applied the MEM technique, which 
is based on Bayes' theorem of probability theory to the 
analysis of the nucleon QCD sum rules. We have investi- 
gated two kinds of sum rules, namely the frequently used 
Borel sum rule and the less known Gaussian sum rule. Be- 
fore analyzing the actual sum rules, we have first tested 
the applicability of the MEM approach by constructing 
and analyzing realistic mock data. Our findings show that 
due to the properties of the kernel and the slow OPE con- 
vergence, it appears to be difficult to extract much mean- 
ingful information on the nucleon ground state from the 
Borel sum rule. Another reason for this failure may also be 
that, because only spectral functions satisfying positivity 
can be analyzed with the currently available MEM proce- 
dure, we cannot use the chiral odd part 772(9^) of Eq.® 
in the present analysis, which has been claimed to be more 
reliable. For instance, in [12] analyses using only 7Ti(g^) 
or 772(9^ ) were carried out and the respective Borel win- 
dows examined. As a result, in case of only using ni{q^), 
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the Borel window is seen to be very narrow, making it dif- 
ficult to obtain a reliable estimate for the nucleon mass. 
On the other hand, in case of using only i72(g^), the borel 
window was shown to be sufficiently large, so that the nu- 
cleon mass can be reliably obtained. It therefore would be 
helpful if n2{q^) could be used. As long as one uses MEM, 
this, however, will only become possible when the parity 
projected sum rules are employed. As an alternative, we 
have formulated the Gaussian sum rule and found that 
it allows us to extract more information on the spectral 
function and enables us to reconstruct the nucleon ground 
state with reasonable precision from both the mock and 
the OPE data. 

As the analysis is done with the MEM technique, we 
obtain the spectral function directly and do not have to 
deal with quantities depending on unphysical parameters 
such as the Borel mass. Moreover, we do not have to re- 
strict the spectral function to the traditional "pole -|- con- 
tinuum" form, which allows us to investigate the spectral 
function of a large variety of of interpolating fields. From 
this investigation, we have found that the nucleon pole is 
independent of the parameter /3 of Eq.Q and vanishes 
when 772 becomes the dominant contribution of the corre- 
lator. Thus we conclude that, in agreement with findings 
of lattice QCD, the nucleon ground state couples only to 
the interpolating field 771 , but not (or only very weakly) to 
rj2- Furthermore, a peak structure is seen around 1.6 GeV 
in the spectral function corresponding to 772, which sug- 
gests that some nucleon resonance in this region couples 
to 772- To clarify the nature of this peak, more thorough 
investigations are needed. 

In all, we have shown that the MEM technique in com- 
bination with the Gaussian sum rule formulated in this pa- 
per is useful for extracting the properties of the nucleon 
ground state and may even make it possible to investigate 
possible excited states. There are, however, still several 
open questions to be answered. First of all, we have so 
far ignored all radiative as corrections to the Wilson co- 
efficients. These corrections are known to be significant 
and it is therefore important to include them for a more 
quantitative analysis. Additionally, possible violations of 
the vacuum saturation approximation should also be con- 
sidered. As a further point, it would be crucial to separate 
the contributions from positive and negative parity states 
to the spectral function, especially to investigate the ex- 
cited nucleon resonances. These issues are left for future 
investigations. 
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